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A novel large-N limit of the multichannel Kondo model is introduced, for representations of the 
impurity spin described by Schwinger bosons. Three cases are found, associated with underscreening, 
overscreening and exact Kondo screening of the impurity. The saddle-point equations derived in this 
limit are reminiscent of the "non-crossing approximation", but preserve the Fermi-liquid nature of 
the model in the exactly screened case. Several physical quantities are computed, both numerically, 
and analytically in the low-oj, T limit, and compared to other approaches. 
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Besides their experimental relevance to magnetic im- 
purities in metals, heavy fermion compounds, and tun- 
neling in mesoscopic systems , quantum impurity mod- 
els provide useful testing grounds for theoretical methods 
dealing with correlated electron systems. Among those, 
explicit solutions in the limit of large spin degeneracy 
(e.g SU(.ZV) spins in the large N limit) have often proven 
to retain many crucial aspects of the low-energy physics 
while being simple to implement 

In this letter, we introduce a novel large-N approach 
to the multichannel Kondo model [||^]. The non- Fermi 
liquid properties of this model in the overscreened regime 
have attracted considerable attention recently. Our ap- 
proach is related to the earlier work of Cox and Ruck- 
enstein ||, in that we also consider the SXJ(N) Kondo 
model with K channels of conduction electrons and take 
the limit N, K — > oo with K/N — 7 being fixed. There 
are two crucial differences with Ref. []5| however: i) we 
are dealing here with Schwinger boson representations of 
the impurity spin, in contrast with the fermionic repre- 
sentations considered in |s| and ii) we keep track of the 
quantum number associated with the "size" of the impu- 
rity spin by imposing a local constraint on the Schwinger 
bosons which is also taken to scale as N. Because of 
these two new features, our approach captures all three 
possible regimes of the multichannel Kondo model (and 
transitions between them) : the underscreened regime (in 
which the Kondo effect only partially screens the im- 
purity spin), the overscreened regime (in which a non- 
Fermi liquid state is formed at low-energy), and an ex- 
actly screened Fermi-liquid regime in between. 

In the large-N limit, our model is solved by a set 
of two coupled integral equations. These equations re- 
semble in structure the "Non-Crossing Approximation" 
(NCA) [|||| , with some crucial differences however asso- 
ciated with points i) and ii) above. Because of our han- 
dling ii) of the constraint, these equations result from a 
true saddle-point principle, with controllable fluctuations 
in l/N. Furthermore, in the exactly screened case, the 
equations derived in this paper preserve the local Fermi- 



liquid character of the problem down to zero-temperature 
and frequency. For this reason, and because Schwinger 
boson mean-field theories yield a quite satisfactory de- 
scription of magnetically ordered phases , our approach 
offers new prospects for a successful treatment of the 
Kondo lattice model. 

We consider a generalised Kondo model with K chan- 
nels of conduction electrons and a spin symmetry group 
extended from SU(2) to SXJ(N). An impurity spin S is 
placed at the origin. We choose to represent the N 2 — 1 
components (S a i3)i< a ,p<N of S in terms of N Schwinger 
bosons b a with a constraint, namely: 

P N 

S a p = bib/3 - — S a(3 , Y h * b * = p C 1 ) 

a=l 

In all the following, the conduction electrons transform 
under the fundamental representation of SU(N). The 
hamiltonian of the model reads: 

ff = E E e P c li* c fn<x + jk s <*p 4i0 c Pia ( 2 ) 

ft Ki<K - a 

P 1< Q <JV pp'iaf3 

In the usual N = 2 case (a, (3 =\,l), this is the mul- 
tichannel Kondo model with an impurity spin of size 
S = P/2. For arbitrary N, Eq. (jl|) means that we have 
restricted ourselves to representations of SU(N) corre- 
sponding to a Young tableau with a single line of P boxes. 
Quantum fluctuations are stronger at small values of P, 
while large P (for fixed N) corresponds to a semiclassical 
limit. 

It is easily checked that a weak antiferromagnetic cou- 
pling (Jk > 0) grows under renormalisation for all K and 
N, and all representations P of the local spin. In order to 
determine whether the renormalization group (RG) flow 
takes Jk all the way to strong coupling (underscreened 
or exactly screened cases), or whether an intermediate 
non- Fermi liquid fixed point exists (overscreened case), 
we have generalised the Nozieres and Blandin stability 
analysis of the strong-coupling fixed point Jk = +00 
Pj. In this limit, a bound-state is formed between the 
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impurity spin and the conduction electrons, which cor- 
responds to a new spin representation dictated by the 
minimisation of the Kondo energy (second term in (0)). 
One must then study the stability of this strong-coupling 
state when a small hopping t of the conduction electrons 
is turned on. The technical steps involved in this anal- 
ysis will be reported elsewhere Kl, and only the main 
conclusions are summarized here. We have established 
that, for arbitrary value of N, there are three possible 
strong-coupling regimes JTo[ ] depending on the number of 
channels K as compared to P : 

• When P > K, (N — 1)K conduction electrons bind 
to the impurity spin in such a way that a free local spin 
of size P sc = P — K is left unscreened for Jjr/i = oo. 
Turning on a hopping tends to increase the total spin of 
the system, corresponding to a weak ferromagnetic resid- 
ual Kondo interaction as one departs from the strong- 
coupling fixed point. Hence, the latter is stable against 
this perturbation, and we have the typical situation of 
an under screened (US) Kondo effect. 

• When P = K, (N — 1)K conduction electrons exactly 
screen the impurity spin and produce a spin singlet state 
(P sc = 0). The strong-coupling state, having the lowest 
possible degeneracy, must be again stable against t/Jx 
and the model displays exact Kondo screening (ES) . 

• When P < K, (N — 1)P conduction electrons screen 
the impurity spin, while (N — 1)(K — P) arrange them- 
selves such that a residual spin P sc = K — P remains. In 
contrast to the above case, the residual Kondo interaction 
produced when a hopping is turned on is now antiferro- 
magnetic, and the strong-coupling fixed point is unstable. 
One expects, as confirmed below, that an intermediate 
fixed point exists with non-Fermi liquid properties: this 
is the overscreened (OS) regime. 

We now turn to the analysis of this model in the large- 
N limit. We shall proceed in a manner similar to Ref. 
||, by setting K = 7 TV and Jk — J/N and taking the 
limit N — > 00 for fixed values of 7 and J. The crucial 
difference (apart from the use of Schwinger bosons) is 
that we shall deal with the constraint in Eq.([j]) by set- 
ting P — pq N and keeping po fixed (instead of fixing 
P = 1 as in ||). By doing so, we are preserving the 
existence of the transition in the large N limit, the three 
regimes above corresponding to po > 7 (US) , po = 7 
(ES) and po < 7 (OS) . This also insures that the model 
is controlled by a true saddle-point at large-N, with con- 
trollable l/N corrections. 

In order to derive the saddle-point equations, we use a 
functional integral formulation of model (|^), with a La- 
grange multiplier field A(r) to implement the constraint. 
Conduction electrons can be integrated out in the bulk, 
keeping only degrees of freedom at the impurity site. The 
local Kondo interaction is decoupled by introducing an 
auxiliary field in each channel Fi(r), conjugate to the 
amplitude J2 a c L( r )^"( T )- This ne ^ wm be responsible 



for capturing the physics of the Kondo effect. Note that 
it is a Grassmanian (anticommuting) field, because of 
our bosonic treatment of the impurity spin. After these 
manipulations, we are left with the effective action: 



,/3 N . p 

S= / dr^bUrJdrbair) + - drJ^F, 
J o a =i J J ° <=1 

+-j= ff\TdT'Y,Fi{r)bUr)G {T-T')F}{T')b a {T') 
+ f^driX{r) hrbi(T)b a (r) - Po Nj 
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In this expression, Go(iu> n ) = J2p — e p) is the on- 

site Green's function associated with the conduction elec- 
tron bath. The quartic term in Eq.(||) can be decoupled 
formally using two bi-local fields Q(t,t') and Q(t,t') 

conjugate to J2 a fra( T )M T ') and J2i F h T ) F i( T ') respec- 
tively. Integrating out all other fields, the action can 
be solved by a saddle-point method over Q, Q and A 
for N — > 00, which leads to coupled equations for the 
Schwinger boson and auxiliary field Green's functions 
G b (r) = - < r&(r)&t(0) >, G f (t) =< TF(r)F t (0) > 
and for the Lagrange multiplier field (the latter is static 
and purely imaginary at the saddle point: iX = A). These 
coupled equations read : 

S 6 (r) = jG (t)G f (t) , E F (r) = G (T)G b (T) (4) 

where the self-energies and Tip are defined by: 



G b V«) = iv n + \- £(,(«/„) 
Gp l {ioj n ) = — - J: F (iu! n ) 



(5) 



In these expressions uj n = (2n + 1)tt//3 and v n = 2mr/(3 
denote fermionic and bosonic Matsubara frequencies. Fi- 
nally, A is determined by the constraint: 



G b (r = 0-) ee J2 G b {iv n )e iv ^ 



Po 



(6) 



We have studied these equations both numerically, and 
analytically in the limit of low-temperature and low- 
energy. The three Kondo regimes are best illustrated by 
Fig.[j], which displays the zero-temperature limit of the 
Curie constant k = limp ->oTximp(T) and of the impu- 
rity entropy S imp = lim T ^ liniy^oo[5'(J 1 ) - S bu i k (T)}, 
as a function of the "size of the spin" po- K vanishes 
in the (OS) and (ES) regimes (po < 7), but reaches a fi- 
nite value K — (po — l)(po — 7 + 1) in the (US) regime 
(Po > 7) 1 corresponding to the Curie constant of a resid- 
ual spin of size P sc = P — K = N(p — 7) in the large-N 
limit. Accordingly, the correlation function of the impu- 
rity spin does not vanish at long times in the (US) regime: 
< S(0)S(t) >~ const., while < S(0)S(r) >~ l/r 2 /( 1+ ^ 
in the (OS) regime. The latter is the behaviour expected 
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from the conformal field theory (CFT) |l3| and Bethe 
Ansatz analysis |l6| of the non-Fermi liquid (OS) fixed- 
point. The residual entropy at T = also takes the 
value expected from the degeneracy associated with a 
free spin of size N(po — 7) in the (US) regime, namely: 
Simp/N = (po~7+!) m (Po - 7 + l)-(Po-7) m bo - 7), 
and vanishes at the (ES) point. In the (OS) regime, its 
value is a universal function of po and 7 associated with 
the non- Fermi liquid fixed point, which will be calculated 
explicitly below (Eq||). 
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FIG. 1. Residual entropy (thin line) and Curie constant 
(bold line) vs. the "size" of the spin po, for 7 = 0.1 (analyti- 
cal expressions given in the text). 
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FIG. 2. Numerically calculated local Curie constant 
Txioc(T) vs. InT with 7 = 0.1, for p = 0.1 ((ES)case) 
andpo = 0.2 ((US) case) 

The cross-over between the high-temperature and low- 
temperature regimes, associated with the Kondo effect, 
is illustrated on Fig which displays the "local Curie 
constant" Txioc(T) vs InT in the (US) and (ES) regime 
(with xioc = < S(0)S{t) > dr). 

We now turn to a more detailed analysis of the low- 
energy, low-temperature behaviour of the above equa- 
tions in the three regimes. In the (OS) case 7 > po, we 



first consider the zero-temperature limit, and perform a 
long-time asymptotic analysis similar to that made in 
Ref. pi. The crucial point in this regime is that the 
constant terms in the r.h.s of Eq.(||) vanish as T — > 0, 
namely: A - T, b (i0+) -> 0,1/ J - Y, F (i0+) ->• 0. As a 
result, a power-law decay G b (r) ~ A b /r 2Ab , Gp(r) ~ 
Ap/t 2Af is found for the Green's functions in the limit 
T^ 1 t (3 — ► 00 (where Tk is the Kondo tem- 
perature). The exponents are given by: 2A& = 1/(1 + 
7) , 2/S.p = 7/(1 + 7), so that the local spin-spin correla- 
tion decays as (S(O)S(r)) cx G b (T)G b {-T) ~ l/r 2 /( 1+ ^. 
The spectral density of the Schwinger boson has the 
following low-frequency behaviour at zero temperature 
(with C a positive constant): 

P^O^^p^ (7) 

It is important to observe that the same power-law 
behaviour is found for u> > and uj < 0, but 
with asymmetric prefactors parametrised by the phases 
6±. The values of these phases can be found using 
a Luttinger-Ward argument leading to: sin# + = 

sin ^(7 -p ) , sin6L = -sin^po- From these 
expressions, we see that Eq.(Q) obeys the positivity 
property appropriate for a bosonic spectral density: 
sign(oj)/9{,(oj) > only as long as po < 7 (P < K). The 
violation of this property signals the transition to the 
(US) regime, in which the solutions just described are no 
longer valid. 

Contact can be made, at least qualitatively, with the 
usual NCA equations (enforcing P = 1 strictly), and with 
Rcf. B, by taking the limit po — > 0, i.e dealing with an 
impurity such that P <C N. In that case, we observe 
that 6- — > 0, and the T = spectral density vanishes 
for negative frequencies. There is no such threshold for 
non-zero values of po however, but the spectral density 
does become increasingly asymmetric as po gets smaller. 

In order to calculate the low-temperature behaviour 
of thermodynamic quantities, the above T = form of 
the spectral densities is insufficient however. The limit 
u)jT — > must be taken while keeping u> = uj/T finite 
10. We have succeeded J|,[)| in calculating analytically 
the spectral functions in this limit, which takes a scal- 
ing form: P6,f(w) = Ab,F T 2Ab - F_1 </Vf(w). The scaling 
function <f>b (found from either a direct solution of the 
integral equations above, or from conformal field theory 
arguments) reads: 



<t>b{u; Po) 



(27T) 



2A h 



2tt 2 



. . S> \T{A b +i 
sinn — ■ 
2 



2tt ' 



(8) 



r(2A b ) 

in which a) = In | sin 6+/ sin^_ |. <j)p has a similar ex- 
pression, with Ap replacing A# and coshw/2 replac- 
ing sinho)/2. Using Eq. (JsJ) and the expression of 
the free energy F imp /N = -p \ - T ^ n hxGb(w n ) + 
7 T En ln ^( lw ") ~ T J2n ^b{iv n )G b (iv n ) we have ob- 
tained [plpl the zero-temperature limit of the entropy, 
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which is a universal number associated with the non- 
Fermi liquid fixed point We find: 



1 + 7 



[/ 7 (l + Po)-/ 7 (l)-/>o)] (9) 



with f-tix) = j 1+ ~' lnsin(u) du. We have also calcu- 



lated pJ9l Si mp for arbitrary finite values of N, K and 
P by applying the conformal-field theory (CFT) meth- 



ods of Ref. Ml3|. We find (for K > P) S r , 



in n; 



7r(jV+ra-l) 



/ sin ■ 



which can also be es- 
p|. The large- 
with cor- 



_j Oiii N + K I 0111 N+K , 

tablished using Bethe Ansatz methods 
N limit of this expression coincides with (f 
rections of order 1/N |l5[ ]. We have also computed the 
low-temperature behaviour of the specific heat ratio and 
impurity susceptibility in the overscreened regime, and 
found : C/T ~ X im P ~ i/T^" 1 )/^) fo r 7 > 1, 
C/T ~ Ximp ~ lnl/T for 7=1, C/T - Xim P ~ const, 
for 7 < 1, also in agreement with the CFT and Bethe 
Ansatz results. Note that the last behaviour holds even 
though the fixed point is not a Fermi liquid (as evidenced 
by the fact that Si mp ^ 0). 

In the (US) regime, the crucial difference with the 
above analysis is that 1/ J — E^(i0 + ) no longer vanishes 
at zero-temperature, but instead diverges logarithmically: 
1/J - E F (iO+) - 1/JelnT + J 2 (while A - E 6 (i0+) ~ 
AiT). This behaviour is well obeyed by the numerical 
solution of the above equations, and has a simple phys- 
ical interpretation: in the (US) case, the residual spin 
is asymptotically free at low temperature, but is weakly 
coupled to the conduction electron gas at any finite tem- 
perature by a ferromagnetic Kondo coupling which van- 
ishes logarithmically as temperature is reduced. Based on 
this observation, we can find the low-temperature scaling 
form of the spectral functions as: Pb(u) = S(u>/Z+\iT) + 
■ ■ ■ and p F (oj) = l/\r?T(j) F {uj/T) + - • •. The form of p b is 
characteristic of an (asymptotically) free local spin. The 
detailed expression of <pF will be given elsewhere ||. 

Finally, we briefly comment on the (ES) case. The 
large-N limit preserves the local impurity properties ex- 



pected from the full Kondo screening, namely «SV, 







and xioc — * const, as T — > 0. The conduction electrons 
form a local Fermi-liquid which is however only weakly 
affected by the screening of the impurity in the large-N 
limit. Indeed, the scattering phase shift obtained from 
Friedel's sum rule is S/tt — 1/N. Accordingly, the scat- 
tering T-matrix, obtained in our approach as the Fourier 
transform of Gf,(r)Gf (— r), is found to have vanishing 
imaginary part in the limit of u>, T — > 0. This is a pe- 
culiarity of the specific spin representation that we have 
considered: considering rectangular Young tableaus with 
P = K = jN columns and N/2 lines would restore max- 
imal unitary scattering S = n/2 for all N. Whether a 
large-N limit exists for that case also is an interesting 



open problem, with potentially useful applications e.g to 
the Kondo lattice model. 

During the course of this study, we learned of simulta- 
neous work by N.Andrei, A.Jerez and G.Zarand on the 
same model using the Bethe Ansatz method fl6|] . Our 
results and conclusions agree when a comparison is pos- 
sible. We are most grateful to N.Andrei for numerous 
discussions. We also thank T.Giamarchi, P.Coleman, 
G.Kotliar, A.Sengupta, A.Ruckenstein and P.Zinn- Justin 
for very helpful suggestions. 
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